Metallic phase of disordered graphene superlattices with long-range correlations 
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Using the transfer matrix method, we study the conductance of the chiral particles through a 
monolayer graphene superlattice with long-range correlated disorder distributed on the potential 
of the barriers. Even though the transmission of the particles through graphene superlattice with 
white noise potentials is suppressed, the transmission is revived in a wide range of angles when the 
potential heights are long-range correlated with a power spectrum S(k) ~ 1/k 13 . As a result, the 
conductance increases with increasing the correlation exponent values gives rise a metallic phase. 
We obtain a phase transition diagram in which a critical correlation exponent depends strongly on 
disorder strength and slightly on the energy of the incident particles. The phase transition, on the 
other hand, appears in all ranges of the energy from propagating to evanescent mode regimes. 



PACS numbers: 73.23.-b,73.63.-b 



I. INTRODUCTION 



The exploration of graphene, a monolayer of carbon 
atoms tightly packed into a honeycomb lattice, has been 
recently attracted special attentions for investigation of 
the fundamental physics and also probable device appli- 
cations such as nanoelectronic based on planar graphene 
structures^* 2 -. In graphene, due to its unique band struc- 
ture with the valence and conduction bands touching at 
two inequivalent Dirac points electrons around the Fermi 
level obey the massless relativistic Dirac equation, which 
results in a linear energy dispersion relation 3 . Massless 
relativistic quasiparticles arising from the cone spectrum 
lead to a number of unusual electronic properties such 
anomalous integer— and fractional^ quantum Hall effects, 
focusing of electron by a rectangular potential barrier 
(Veselago lensing)^, special Andreev reflection 7 -, observa- 
tion of the plasmaron composite- and minimal conduc- 
tivity 9 . 

Interestingly, relativistic quantum quasiparticles inci- 
dent normally to a high electrostatic potential barrier in 
graphene can pass through it with perfect transmission 
regardless of the height and width of the barrier—. This 
phenomenon which is referred as the Klein tunneling is 
in contrast with the quantum massive carrier tunneling 
where the transmission probability decays exponentially 
with increasing of the barrier height and width. Recently, 
evidences for the Klein tunneling of the Dirac fermions 
across p-n junction have been experimentally observed 
when a gate-induced potential step is steep enoug h 11 ' 12 . 

In graphene sheets the type of particle (electrons or 
holes) and the density of the carriers can be controlled 
by tuning a gate bias voltage^ 3 . - — . Moreover, graphene 
superlattices may be fabricated by adsorbing adatoms on 
the graphene surface by positioning and aligning impuri- 
ties with scanning tunneling microscopy^, or by apply- 
ing a local top gate voltage to graphene^ 7 -. The transition 
of hitting massless particles in a cleani^ or disordered^ 
graphene-based superlattice structure has been studied. 
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FIG. 1: Graphene superlattice with long-range correlated dis- 
order on the potential barriers. 



It is shown that the conductivity of the system depends 
on the superlattice structural parameters. 

The first study on electronic properties of monolayer 
and bilayer graphene superlattices was performed by Bai 
and ZhangAS. They showed that the angularly averaged 
conductivity can be controlled by changing the struc- 
ture parameters. It has been shown that massless Dirac 
fermions are generated in one-dimensional external peri- 
odic potential close to the original Dirac poin t 20 ' 21 . The 
Dirac points depend on geometrical parameters for in- 
stance the potential of the barriers/wells, the period of 
the potential and transverse wave number—. An evidence 
for such Dirac points is the conductance resonances are 
appeared in the special potential values 2 ^. Moreover, the 
conductance of graphene superlattice with uncorrelated 
disorder on the width of the barriers was calculated in 
[ [l9| . It was shown that the transmission of the quasi- 
particles with large angles incidence to the potential bar- 
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riers is suppressed by disorder strength and the sample 
size too. Therefore, the results of the finite-size scaling 
computations predicted a zero conductance for all the 
graphene super lattices, except for some resonant barrier 
thickness in which the conductance tends to a nonzero 
constant in the thermodynamic limi^. 

A number of numerical calculations of electron trans- 
port confirmed the absence of the localization in the pres- 
ence of the long-range random potential in disordered 
graphene^. The main quantity mostly studied numeri- 
cally is the conductance, G of a finite-size graphene sam- 
ple with a width, W much larger than the length, L. The 
setup allows us to define the "conductivity" a — GL/W 
even for ballistic samples with L much shorter than the 
mean free path, I. 

It is well known that the transmission of the quantum 
massive carriers unexpectedly increases when special cor- 
relation is applied on disorder^. This is in contrast with 
the Anderson localization in which all states are expo- 
nentially localized in one dimensional uncorrelated dis- 
order. Experimental evidence for discrete the number of 
the extended states has been observed in random-dimer 
semiconductor superlattices as a short-range correlated 
disorder 25 . Long-range correlated sequences of the po- 
tential barriers in semiconductor superlattices, however, 
could result in a continuum of extended states giving rise 
the mobility edges 2 ^— . 

In this paper, we study the conductance of massless 
Dirac fermions through graphene superlattices with a 
long-range correlated disorder on the potentials of the 
barriers. Transmission of the large angles incident elec- 
trons to graphene superlattice which are suppressed by 
uncorrelated disorder, are revived by applying correlation 
between random potentials of the barriers. As a result, 
the conductance increases with the correlation exponent. 
Consequently, an insulator-to-metal transition emerges 
at a critical correlation exponent which depends strongly 
on the disorder strength. One should notice that such 
phase transition emerges for all ranges of the energy. In 
addition, the dependence of the conductance to the su- 
perlattice parameters is investigated for different corre- 
lation strengths. 

The paper is organized as follows: we present the trans- 
fer matrix and Fourier Filtering method used to calcu- 
late transmission through long-range correlated graphene 
superlattices in Section II. Our results and discussions 
will be presented in Section IV presenting the metal-to- 
insulator phase transition and investigating the emer- 
gence of the phase transition along different energy 
ranges. The last Section concludes our results. 



II. MODEL AND THE TRANSFER MATRIX 

In the low-energy limit, charge carriers near the Dirac 
point in the continuum model can be described by the 
following non-interacting Hamiltonian 
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FIG. 2: Random distribution of the correlated sequences gen- 
erated by Fourier Filtering Method, a) uncorrelated case cor- 
responding to q = 0.500, b) a = 1.766, c) a = 1.993. 



H = -ihv F ~^.^ + V(x) (1) 

where Vp — 10 6 m/s is the Fermi velocity and a = 
(cr X) (T y ), are Pauli matrices. We consider a lattice of the 
electrostatic potentials as barriers which are induced by 
top gate voltages. Therefore, the potential of the barriers 
are sorted as the following: 

!Vi d%i-\ < x < d 2i , i = l,2, ... 
(2) 
elsewhere 

A schematic representation of graphene superlattice is 
shown in Fig. [I] We consider the number of regions (wells 
and barriers) to be equal N and then there is (N — l)/2 
barriers. The width of the barriers and wells are consid- 
ered to be fixed. The height of the barriers fluctuates 
around its mean is defined by Vi — (V)(l + <r £j) where a 
is the variance of the potentials of the barriers and {e^} 
is a long-range correlated random sequence of data with 
the Gaussian distribution. The random sequences will 
become normalized, accordingly the mean value £j is set 
to be zero and its variance has been fixed. 

Before passing to calculate the conductance, we might 
generate a random sequence with long-range correlation, 
{ei} will be considered to describe the trace of a fractional 
Brownian motion with a power spectrum S(k) ~ \/k@ 
where 1 < < 3 and = 2H + 1 = 2a - 1 (1 < 
a < 2). Here H is the Hurst exponent. In the case 
of power-law decaying auto-correlations, the correlation 
function decays with an exponent 7 such that C(xi — 
Xj) oc| i — j |~ 7 where 7 = 1—0. Random sequences with 
weaker positive correlation is generated with 0.5 < a < 1 
which is referred to fractional Gaussian noised. In this 
case, = 2H — 1 = 2a — 1 and a — 0.5 corresponds to 
uncorrelated disorder or white-noise. 
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FIG. 3: Transmission in terms of the incident angle hitting 
to graphene superlattice with various correlation exponents 
a. The number of the barriers is N = 1000. Averaged po- 
tential of the barriers and energy of the incident electrons are 
(V) = 200meV and E — 50meV, respectively. Here, since 
£ = —3, and transmission is calculated in the presence of 
purely propagating modes. 



It is commonly usual to apply Fourier filtering 
method 2 ^ to generate a sequence with the long-range 
correlation. The method is based on a transformation 
of the Fourier components (6k) of a random number se- 
quence {0i}. Uncorrelated random numbers {6i} have a 
Gaussian distribution. Finally, the inverse Fourier trans- 
formation of the sequence {£&}(= k~^9k) results to the 
interested sequence {e,}. Three landscapes of random 
data generated by the mentioned method are shown in 
Fig. [2] for different correlation exponents. Clearly, corre- 
lation between random potentials leads to a reduction in 
the fluctuations of the random distribution. We should 
notice that we have checked the invariance of our results 
in comparison with those results obtained from random 
sequences produced by the mid-point method 30 . 

Having such configuration of the potentials gives rise a 
superlattice consists of two types of graphene as electron- 
doped or hole-doped. The doping type of graphene in 
wells is n-type for quasiparticles with incident energy 
E = hvp/X > where A is electron wavelength in wells. 
The doping type in barriers depends on the potential 
height Vi- In other words, the i th barrier is hole-doped 
graphene if E < Vi, while it is electron-doped if E > Vi. 
Accordingly, the type of doping in barriers is also altered 
randomly. We assume that the angle of electrons inci- 
dence to superlattice is <f> — ip± along the x axes. How- 
ever, the transfer matrix takes the angle of each region 
into account separately. The general solution of Eq. ([1]) 
results in the following spinor for the i th region; 



ip(x,y) = a. t 



i{k i:c x+k y y) 



Sie 



i(-TT — ipi) 



(3) 



where and bi are the transmission and reflection am- 
plitudes, respectively. Other parameters in the spinor 
are 



Si = sgn(E - Vi), k y — kp sin(tpi) = k F sin0 
<J((E -V{x))/hv F ) 2 -k 2 y , Vi = arctan(Vfc«)-( 4 ) 



If the energy of the incident electrons being close to 
Vi value in the i th barrier, ki X will become imaginary 
value in some angles, resulting in an evanescent mode. 
In disordered graphene superlattice, evanescent modes 
emerge when E ~ (V). The transfer matrix is extracted 
by the continuity of the wave functions at the junction 
interfaces. It can make a relation between wave functions 
of two sides of a step potential from the i th region to the 
(i + l) th like 



bi+i 



= Mi. 



where the transfer matrix Mj_i_i , is; 



(5) 



Mi+i, 



ran mi2 
m* 12 m n 



and the matrix elements of M are as follows; 



v Zs%+\ cos tfi+1 

ra 12 = e -i(fc»+fc(i+i).)^( a '+i e ~' yi+1 -* e ~ iyi ) 

V 2si+i cos (pi+t * 



(6) 



) (7) 



The current in the x— direction and in the i th region 
can be derived as the following; 



J l x = v F ^a x ^ = 2v F Si cos ifi(\ at | 2 - | bi | 2 ) (8) 

Current conservation between regions i th and j th im- 
plies on 



I «i I 2 - I b o | 2 = Det[M jti ]{\ ai | 2 - | bi | 2 ) (9) 



where 



Det[Mj t i\ = Si cos Lfi I Sj cos <pj 



(10) 
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FIG. 4: Transmission as a function of the incident angle 
through graphene superlattice with long-range correlated ran- 
dom potentials for various barrier numbers N. Here a = 1.102 
and a = 0.1. 



FIG. 5: Conductance through graphene superlattice with 
long-range correlated random potentials as a function of bar- 
rier number N for a = 0.1. This is a finite size scaling for 
different correlation exponents a. Here, (V) = 200meV and 
E = 50meV. 



The total transfer matrix which makes a relation be- 
tween incident and transmitted wave functions is a se- 
ries product of the transfer matrices arising from each 
interface. For N regions incorporating of the barriers 
and wells, regarding to Eg. dill) , the matrix is denned as 
P = M N>1 . 

If the first and last regions of superlattice are electron- 
doped graphene, transmission probability for (N — l)/2 
barriers can be calculated by means of the product ma- 
trices for 6/v = as the following 



T (E t 4>) = ^ = —^ — , (11) 
1 ,w J} n Det[P] oi ' { ' 

where J^ ut and J} n are out- and inflowing currents, re- 
spectively. Because the configuration of the potential 
barriers is considered such a way that s\ = sn and 
"Pi = <Pn = 0j the- conservation of the current be- 
tween the first and last regions implies that Det[P] = 1. 
Therefore, the transmission formula can be simplified as 
T{E,4>) = 1/ | P22 | 2 where a N / ai = Det[P]/P 22 . Fi- 
nally, using Landauer-Biittiker formula^! and an angu- 
larly averaging, the conductance is obtained by the fol- 
lowing integration. 



,tt/2 

G = G T(E,<j>)co&{4>)d(t> (12) 

J -it/2 

where Go = e 2 mv-pW/h 2 . Here W is the finite width of 
graphene ribbon along the y— direction. 



III. RESULTS AND DISCUSSION 

A. Phase transition 

Let us first calculate the transmission probability and 
study the electronic properties of disordered graphene 
superlattices. The transmission of electrons hitting to 
a disordered graphene superlattice as a function of the 
incident angle is shown in Fig. [3] for several values of 
correlation strengths characterizing with the correlation 
exponent a. In all calculations, barrier and well widths 
are considered to be D = 50nm and L — 30nm, re- 
spectively. Moreover, we assumed that the energy of the 
charge carriers and the averaged potential of barriers be- 
ing E = 50meV and (V) = 200meV, respectively. The 
wavelength of the incident electrons is thus A = 83nm. 
Therefore, with such parameters we surely conclude that 
the transmission of the charge carriers shown in Fig. [3] 
(with £ = (E — (V))/E — —3) is a purely propagat- 
ing mode. As shown in this figure, the transmission of 
the electrons hitting to superlattice with large angles in- 
creases by increasing correlation between the random po- 
tentials of the barriers. In other words, applying corre- 
lation between random potentials of the barriers causes 
to extend the angular window of the conducting mode 
around the normal incidence. This effect is in contrast 
with those results obtained in uncorrelated potentials of 
the barriers in which the transmission of the massless 
carriers is suppressed for the large incident angles except 
at <f) — 0. Perfect transmission at normal incidence can 
be described by the Klein tunneling. Similarly, in the 
correlated case, by increasing disorder strength and also 
the number of the barriers, the transmission of the quasi- 
particles is suppressed at all ranges of the incident angle 
except at <fi = 0. Fig. 0] shows the suppression of the 
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FIG. 6: Exponent of the conductance in a power law form 
G/Gq oc N~ v as a function of correlation exponent a. In 
the inset the fitting of a Fermi-Dirac like function (Eq |14|) to 
numerical data is shown. 
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FIG. 7: Metal-to-insulator phase diagram. The critical cor- 
relation exponent, a cr increases with disorder strength, a. 



transmission at large incident angles when the number 
of the potential barriers increases. 

In linear regime, the conductance is proportional to an- 
gularly averaged transmission projected along the current 
direction. To understand how the correlation between 
random potentials affects on the transport properties of 
graphene superlattice, we calculated the size dependence 
of the conductance for the various values of the corre- 
lation exponents and results are plotted in Fig. [S] It is 
clear from the figure that there is a critical correlation 
exponent value of a such that for a < a cr , the con- 
ductance decreases with increasing system size, while it 
goes to a constant value for a > a cr . In other words, 
a consequence of applying long-range correlation is the 
emergence of a phase transition from the insulating to 
metallic states. It is worthwhile noting that this phase 
transition is not a finite size effect. The conductance de- 
creases with the number of the barriers as a power law 
behavior. The following function is fitted to a log-log 
plot of the conductance. 

— oc N-^ a '^ (13) 
Go 

where rj as functions of a and a decreases by increasing 
the correlation exponent. Furthermore, disorder strength 
strongly suppresses the conductance such that the critical 
correlation exponent which implies on the emergence of 
a phase transition, increases with the disorder strength. 

Conductance through disordered potentials of barri- 
ers having correlation exponent a = 0.5 decays with 
the barriers number as iV~ 0,5 . The same decaying of 
the conductance was reported in Ref. [l9| for graphene 
superlattice with white noise potentials distributed on 
the barriers. However, applying a long-range correlation 
between random potentials facilitates the conductance 



through graphene superlattice. 

A power-law fitting for the conductance is shown in 
Fig. [SJ It is determined a dependency of rj on the corre- 
lation exponent a and disorder strength a which is rep- 
resented in Fig. [5] For the sake of having the critical cor- 
relation exponent function, we have used a fit function 
of the exponent function rj(a, a) in a transition function 
such as a Fermi-Dirac function: 

[ a » a cr 



(14) 

where /3 and a cr are fitted parameters and 7 = 0.5 for £ = 
—3. The range that 77(a) decreases from the value of 7 to 
zero is related with the inverse of /3. The emergence of a 
transition from insulating to metallic phase corresponds 
to the variation of 77 from 7 to 0. Now, we provide a 
phase transition diagram in which the critical correlation 
exponent depends on disorder strength. Fig.[7]shows that 
the critical correlation exponent increases when disorder 
strength increases up to a — 0.3. Roughly speaking, the 
system remains in a metallic phase when the correlation 
is very long range at low disorder strength whereas it 
turns out to be an insulator at large disorder strength 
values and at mid-range correlation. 

B. Energy range of the phase transition 

To provide a fascinating experimental manifestation of 
the phase transition, it is significant to demonstrate that 
the phase transition can exist in a continuum range of 
energies not just at some discrete energies. 

Let us now concentrate on the conductance behavior 
in the different range of the Fermi energy. For a sin- 
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FIG. 8: Functional of the conductance in terms of £ which 
represents Fermi energy. There are two peaks in conductance 
correspond to resonant states. In the inset 3D contour plot 
of the transmission in the plane of £ and the incident angle. 
Resonant states corresponds to more angularly open a domain 
for transport. Here N = 800 and a = 0.1. 




FIG. 9: a) Conductance in terms of the barrier numbers for 
graphene superlattice with white-noise disorder on potentials 
of the barriers and for different values of the Fermi level £ . b) 
Exponent of the conductance (rj) in a power-law form G/Go oc 
N~ v as a function of £. 



gle barrier on graphene, in the range of —1 < £ < 1, 
it is proved that both the evanescent and propagating 
modes coexist^. Out of this range, all states are fully 
in the propagating modes. By considering this fact, we 
investigate the conductance as a function of £ through 
a disordered graphene superlattice for several correlation 
exponent values and the results are shown in Fig. [5J It 
can be seen that the same behavior as a single barrier 
case is appeared in the different ranges of £. Contri- 
bution of the evanescent modes in the conductance is 
dominant in £ = and therefore the conductance is sup- 
pressed in this point. For £| — > 1, the contribution of 
the propagating modes becomes dominant and thus the 
conductance increases. For £ < 0, the conductance oscil- 
lates at resonant states which originates from the perfect 
tunneling of the charge carriers near to a normal inci- 
dence. A 3D contour plot of the transmission in terms 
of £ and the incident angle (see the inset Fig. [5J shows 
that at resonant states, a conducting domain of the an- 
gles is opened around the normal incidence. The reso- 
nant condition for a single barrier in graphene is given 
by k-pDyJ^ 2 — u 2 = mr where u = sin ip. Distance of 
the resonant states can be extracted from expanding the 
resonant condition for normal incidence u <C £. In this 
case, two sequential resonant states have a distance like 
A£ r = ^7d(1 — le?)- As a result, the period of the con- 
ductance reduces with D, L while enhances with |£|. In 
the limit of £ —> 1 which it means £ > < V >, the 
conducting channels are opened for the whole range of 
the angles, T(<j>) = 1, and thus by using Eq. ([T2]). we 
then get lirri£_>i G/Go = 2. 

By applying a long-range correlation between the po- 
tentials of the barriers, Fig. [8] indicates that the con- 
ductance shows an enhancement in some ranges of the 
energy compared to the resonant states and the range 



including the evanescent modes | £ |< 1. For clarifying 
that, firstly we investigate the conductance suppression 
with the number of the barriers over all ranges of the 
Fermi energy for a sequence of white noise disorder dis- 
tributed on the potentials of the barriers. Results are 
represented in Fig. which show a power-law form of 
the conductance in terms of the barrier numbers for dif- 
ferent values of £. It is significant to understand how the 
exponent of r\ varies with the Fermi energy. The slope 
of the lines in log-log plot of Fig. [5^ for different £ are 
indicated in Fig. \§]p. It is clear that close to resonant 
states £ r = —1.85, —0.45, the exponent of r\ reduces from 
the value of 0.5 to its resonant values r\ r — 0.12,0.25, 
respectively. As a consequence, close to resonant states, 
the suppression of the conductance which is induced by 
random potentials of the barriers is much weaker than 
other states. The resonant value of r\ r decreases when 
£ goes away from the region including the evanescent 
modes | £ |< 1. 

Our numerical calculations demonstrate that the 
insulator-to-metal phase transition occurs in all ranges 
of the energy. It is clear from Fig. [TU] that the phase 
transition appears not only at the propagating and reso- 
nant states, but also in the fully evanescent mode £ = 0. 
In Fig. [TQl Tj{a — 0.5) decreases to 0.3 at £ = —2 which 
is close to the second resonant state shown in Fig. [5J An- 
other exception around the resonant states is the width 
of function 77(a) which increases at the resonant states. 
Therefore, at the resonant states, the transition from in- 
sulating to metallic phase is smooth along the correlation 
exponent. The inset Fig. [TQl shows small fluctuations of 
the critical correlation exponent as a function of £. Ac- 
cordingly, the phase transition is a general behavior of 
all ranges of the energy. 
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FIG. 10: Exponent of the conductance in a power law form 
G/Go oc N~ v as a function of the correlation exponent a for 
different values of f . In the inset the metal-to-insulator phase 
transition in all ranges of the energy. 

C. Resonance in the conductance 

Now, we study the effect of the long-range correlated 
disorder on the resonance phenomena seen in the conduc- 
tance. Resonance in the conductance of graphene super- 
lattice with white noise disorder distributed on width of 
the barriers has been studied before^. Fig. [Til shows the 
conductance oscillations as a function of the barrier width 
(D) and also distance between barriers (L) for several 
values of the correlation exponent. Apparently, the ap- 
plication of the long range correlation between potentials 
of the barriers increases the conductance for all ranges of 
D and L. Moreover, as shown in Fig. ITTa . the conduc- 
tance tends to a constant value after some oscillations in 
thin barriers independent of the barrier width. In fact, 
the transmission of the quasiparticles hitting to graphene 
superlattice at large incident angles is suppressed for the 
wide barriers, and thus only transmission arising from the 
Klein tunneling around the normal incident contributes 
into the conductance integration given by Eq. (TT21) . The 
same behavior occurs in the case of N = 1. Configu- 
ration average of the transmission through one barrier- 
depends on the width by functions of < sm 2 (k x D) >c.A. 
and < cos 2 (k x D) >c.A. where k x is a random parame- 
ter. It is trivial that the average of the transmission and 
consequently the conductance become independent of the 
width when D value increases. On the other hand, by ap- 
plying of the randomness on k x , the resonant condition 
k x D = mr in wide barriers can not be satisfied, and there 



is thus no longer the resonant peak in the conductance. 

Fig. ITTb shows the conductance oscillations as a func- 
tion of L. In this case, the decaying of the oscillations 
is much weaker than the conductance oscillations with 
the barrier width. In fact, since there is no disorder in 
the wells, in a fixed barrier width, the resonant condition 
affects much less than the barrier resonant condition. 
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FIG. 11: Conductance oscillations in terms of a) the width 
(D) of the barriers b) the distance (L) between barriers for 
different correlation exponents. Here N = 800 and a — 0.1. 



IV. CONCLUSION 

By using the transfer matrix method, we investigate 
the conductance through a graphene superlattice with 
long-range correlated disorder distributed on potentials 
of the barriers. Applying a correlation between poten- 
tials opens angularly domain window of the conducting 
channels in competition with the factor of the disorder 
strength which suppress the transmission at large inci- 
dent angles. As a result, the conductance increases with 
the correlation between the potentials of the barriers 
gives rise a metallic phase. We obtain a phase transi- 
tion diagram in which the critical correlation exponent 
for such a phase transition depends strongly on disor- 
der strength and slightly on the energy of the incident 
particles. At resonant states, the suppression of the con- 
ductance with the number of the barriers is much less 
than other states. Our finding for the dc conductance 
of the graphene superlattices should be important to the 
design of electronic nanodevices based on graphene su- 
perlattices. 
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